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ABSTRACT 

Stars form in dense cores of molecular clouds that are observed to be signifi- 
cantly magnetized. A dynamically important magnetic field presents a significant 
obstacle to the formation of protostellar disks. Recent studies have shown that 
magnetic braking is strong enough to suppress the formation of rotationally sup- 
ported disks in the ideal MHD limit. Whether non-ideal MHD effects can enable 
disk formation remains unsettled. We carry out a first study on how disk forma- 
tion in magnetic clouds is modified by the Hall effect, the least explored of the 
three non-ideal MHD effects in star formation (the other two being ambipolar dif- 
fusion and Ohmic dissipation). For illustrative purposes, we consider a simplified 
problem of a non-self-gravitating, magnetized envelope collapsing onto a central 
protostar of fixed mass. We find that the Hall effect can spin up the inner part of 
the collapsing flow to Keplerian speed, producing a rotationally supported disk. 
The disk is generated through a Hall-induced magnetic torque. Disk formation 
occurs even when the envelope is initially non-rotating, provided that the Hall 
coefficient is large enough. When the magnetic field orientation is flipped, the 
direction of disk rotation is reversed as well. The implication is that the Hall 
effect can in principle produce both regularly rotating and counter-rotating disks 
around protostars. We conclude that the Hall effect is an important factor to con- 
sider in studying the angular momentum evolution of magnetized star formation 
in general and disk formation in particular. 

Subject headings: accretion, accretion disks — magnetic fields — ISM: clouds — 
stars: formation 



1. Introduction 

Disks are an integral part of star formation; they are the birthplace of planets. How 
they form is a long-standing, unresolved problem. A major difficulty is that their formation 
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is greatly affected by magnetic braking, which has been hard to quantify until recently. 



There is now increasing theoretical evidence that magnetic braking may suppress the 
formation of rotationally supported disks (RSDs hereafter) in dense core s magnetized to a re- 
alisti c l evel, with di r nensio nless mass-to-fiux ratios A of a few to several (ITroland fc Crutcher 
20081 ). lAUen et al.l (120031 ) first demonstrated through 2D (axisymmetric) simula tions that 
RSDs are suppressed by a moderately strong magnetic field in the ideal MHD limit. iGalli et al 
(120061 ) showed analytically that the disk suppression is due to the formation of a split mag- 
netic monopole, which is an unavoi dable consequence of fl ux freezing. The efficient disk 



braking was co nfirmed numerically b v lMellon fc Lil (120081) and 



see, however, iMachida et al.l |2010| and iHennebelle fc Ciardi 



Hennebelle fc Fromand (120081 ) 
2009 for a different view, and 



Li et al.ll201ll for a more detailed discussion). The ideal MHD approximation must break 
down in order for RSD to form. 

Two non-ideal MHD effects have alread y been explored in the context o f disk formation: 
ambipolar diffusion and Ohmic dissipation. iKrasnopolsky fc Konigll (120021 ) found that am- 
bipolar diffusion tends to make disk braking more efficient, because it enables the magnetic 
flux that would have gone into the central protostellar object (and form a split monopole) in 
the ideal MHD limit t o pile up at small radii outside of the object instead, making the regio n 
strongly magnetize d (ILi k. McKed Il996l . ICiolek fc Konigll Il998l . IContopoulos et al.l Il998l ) . 
Mellon fc Lil ( l2009l ) showed that the enhanced braking is strong enough to suppress the for- 
mation of RSDs for realistic core condition s. The effect of Ohmic dissipation was examined 
by several groups, starting with IShu et al.l (120061 ). They suggested that enhanced resistiv- 
ity (well above the classical value) is needed for Ohmic dissipation to weaken the magnetic 
braking enou gh to form a relative l y larg e RSD of tens of AUs or mor e. The suggesti o n was 
confirmed by iKrasnopolsky et al.l (120101 : KLSIO hereafter), although iMachida et al.l (120101) 
were apparently able to form RSDs u sing the classical resistivity computed by iNakano et al. 
(120021 ) (see also iDapp fc Basull2010l ). In any case, a third non-ideal MHD effect, the Hall 
effect, h as not been e xplored in the context of disk formation (see, however, the independent 
work of lC. Braiding in her unpublished PhD thesis); it is the focus of this paper. 



Our goal is to determine whether the Hall effect by itself can enable an RSD to form 
in the presence of a relatively strong magnetic field and, if yes, to estimate the magnitude 
of Hall coefficient needed for RSD formation. We find that the Hall effect can actively spin 
up the inner part of the protostellar collapsing flow, potentially to Keplerian speed, unlike 
the other two non-ideal MHD effects. The combined effect of all these three non-ideal MHD 



terms on disk formation will be explored in another investigation (ILi et al.l 1201 ll : LKSll 
hereafter) . 
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2. Problem Setup 



We adopt the same problem setup as in KLSIO, where we consider the collapse onto a 
star of 0.5 M0 of a rotating envelope that is uniformly magnetized initially. In the absence of 
any non-ideal MHD effect, the magnetic field would prevent an RSD from forming through 
magnetic braking. The Hall effect changes the evolution of the magnetic field (and thus the 
braking efficiency) through the induction equation: 

OB 

'dt 



V(v X B) - V X {g[(V X B) X B]} 



T. 



where Q is a coefficient of the Hall effect which, for simplicity, we will assume to be spatially 
constant . The MHD equations are solved using ZeusTW, a 3D non-ideal MHD code based on 



Zeus3D fjClarke. Norman. Sz Fiedler Il994r). We tre a t the Hall t erm in the in duction equation 



Sano fc Ston3 J2002h and [Hubal J2003h . which includes 



using an explicit method based on 
subcycling to speed up the computation. 

As in KLSIO, we adopt a spherical polar coordinate system (r, 6, 0), and fill the 
computation domain with an isothermal gas of sound speed a = 0.2kms~^ between rj = 
1.5 X 10^^ cm and = 1.5 x 10^^ cm. A uniform density po = 1-4 x 10~^^gcm~^ is assumed, 
so that the total envelope mass is IMq. For the initial rotation, we adopt the following 
prescription: 

v<p = U0,o tanh(cc7/ roj (2) 

where w is the cylindrical radius, and v^^o = 2 x lO^cms"^ and Wc = 3 x 10^^ cm are 
chosen. With this setup, an RSD of ~ 400 AU in radius is produced at a representative time 
of t = 10^^ s in the absence of any magnetic field (see Fig. 1 of KLSIO). It is completely 
destroyed by a moderately strong initial magnetic field of Bq = 35.4 /xG in the ideal MHD 
limit (see Fig. 2 of KLSIO). Whether the RSD can be restored by the Hall effect is the 
question that we seek to address. 



3. Result 

We consider first a model where the initial magnetic field is parallel to the rotation axis 
(Model PARA in Table [1]) and the coefficient Q' = 3 x 10^^ in the Lorentz-Heaviside units 
convenient to use with the Zeus family of codes, corresponding to Q = Q' /y/An = 8.46 x 10^^ 
in Gaussian CGS units (cm^s~^G~^). An RSD is able to form in this illustrated 
Fig. [H which plots a color map of density distribution, velocity field and magnetic field lines 
in the meridian plane (left panel), and the infall and rotation speeds on the equator (right 
panel), at a representative time t = 2 x 10^^ s. The RSD shows up clearly on the color map 
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as the flattened high density region of about 60 AU in radius. It has an infall speed that 
is close to zero and a rotation speed that is Keplerian. The Hall effect has clearly enabled 
an RSD to form, although it is not by weakening the magnetic braking that would have 
suppressed the RSD formation in the ideal MHD limit, unlike the case of Ohmic dissipation 
(KLSIO). The RSD is formed because the inner part of the collapsing flow is actively spun 
up by the Hall effect. A strong support for Hall spin- up as the cause for the RSD formation 
comes from a second simulation where we flip the sign of the initial magnetic field direction 
(Model ANTI), so that the field is anti-parallel, rather than parallel, to the rotation axis. In 
this case, an RSD disk is again formed, but it rotates in a direction opposite to the initial 
rotation, as shown in Fig. |2l The reversal of the disk rotation direction means that the Hall 
effect does not merely weaken the magnetic braking so that enough of the original angular 
momentum is retained for RSD formation. 

The reason that the disk rotation direction is reve rsed when the magnetic field direction 



is fiipped is the following (see lWardle &: Konigllll993l for a similar discussion in the context 



of disk- wind launching). From the induction equation ([T]), it is easy to see that the Hall 
effect forces the field lines to move with a velocity 

VH = -Q(V X B) = -^J (3) 

c 

where J is the current density. The pinching of the poloidal magnetic field near the equatorial 
region produces a strong toroidal current, j^, which changes sign as the poloidal field direction 
is fiipped. The toroidal current of different signs causes the field lines to wind up toroidally 
in opposite directions, giving rise to magnetic torques of opposite sense. In the case where 
the initial field line is parallel to the rotation axis, is positive, giving rise to a toroidal field- 
line speed fH,,/, that is ne gative (assuming th at the Hall coefficient is positive, which is not 



necessarily true, see, e.g.. lWardle fc Nglll999l and LKSll; wh,0 changes sign if Q is negative). 
Since the toroidal current is the strongest on the equator, the field line is twisted backward 
most on the equator and less so above and below the equatorial plane. The differential twist 
forces the field line to bend in the positive azimuthal direction, creating a torque that spins 
up the material near the equatorial plane in the same direction as the initial rotation. This 
explains why the RSD in the field-rotation aligned case rotates in the same sense as the 
initial rotation. When the initial field direction is fiipped, the toroidal current (and thus 
V}i,ip) changes sign, forcing the field lines to bend in the negative azimuthal direction. The 
resulting negative torque not only overcame the initial positive rotation but also produced 
a counter-rotating RSD in the anti-aligned case. 

The Hall spin-up is even more unambiguous in Model NoROT where the envelope is 
initially non- rotating; any subsequent rotation must be due to the Hall effect. The right panel 
of Fig. [3] shows that an RSD is again formed in this case, with the rotation speed approaching 
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loQioP (cgs units) 
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Fig. 1. — Left: Snapshot of the rotationally supported disk formed through Hall effect in 
the meridian plane at a representative time t = 2 x 10^^ s for Model PARA. Logarithm of 
density (color map); radial infall sonic transition Vr = —a = —2 x lO^cms"^ (white line); 
poloidal velocity field (arrows); magnetic field lines (gray solid lines); level of Keplerian 
support \v^/vk\ (dark gray dashes at 90%), where vk = {GM^:/wY^'^. Right: infall and 
rotation speeds on the equator. Also plotted for comparison are the Keplerian speed and 
free fall speed, the sound speed (horizontal dashed line), and zero speed line. 
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Table 1. Model Parameters 



Model 


Q [cgs units) 


Bo ifiG) 


Rotation 


RSD? 


PARA 


8.46 X 10^1 


35.4 


yes 


yes 


ANTI 


8.46 X lO^i 


- 35.4 


yes 


yes 


NoROT 


8.46 X lO^i 


35.4 


no 


yes 


LoB 


8.46 X lO^i 


10.6 


no 


yes 


LoQ 


8.46 X 10^0 


35.4 


no 


yes 


LoQl 


2.82 X 10^° 


35.4 


no 


yes 


LoQ2 


8.46 X 10^9 


35.4 


no 


no 




Fig. 2. — Left: Snapshot of the rotationally supported disk formed through Hall effect in 
the equatorial plane at a representative time t = 2 x 10^^ s in the equatorial plane, with 
anti-aligned initial magnetic field and rotation axis (Model ANTI). Note that the Hall effect 
forces the inner part rotate in a direction opposite to that of the outer part. Right: infall 
and rotation speeds on the equator, with the Keplerian and free fall speeds also plotted for 
comparison. 
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the Keplerian speed and the infall speed dropping close to zero inside a radius ~ 10^^ cm. 
Since the envelope is initially non-rotating, the spin of the RSD must be balanced by material 
rotating in the opposite direction. The left panel of Fig. [3] demonstrates that this is indeed the 
case. It shows that, outside the equatorial plane, materials of positive and negative angular 
momentum occupy alternating shells. Inside a radius of ~ 5 x 10^ AU, the angular momentum 
is positive, dominated by the (small) RSD and highly flattened pseudodisk outside RSD (see 
the isodensity contours). At larger distances, the angular momentum becomes negative, 
dominated by the outermost counter-rotating shell shown in the Fig. [31 The total positive 
angular momentum (1.73 x lO^^gcm^s^^) does not exactly balance out the total negative 
angular momentum (—2.12 x 10^^ gcm^ s~^), however, indicating some angular momentum 
is transported out of the simulation box, presumably through torsional Alfven waves. Note 
that the Hall effect modifies the envelope rotation well beyond the 10^ AU-scale RSD. We 
have verified that when the initial magnetic field direction is flipped, the direction of the 
Hall-induced rotation is reversed. 

Besides torquing up the envelope in the toroidal direction, the Hall effect also enables 
the magnetic field lines to diffuse out radially. In the absence of any radial diffusion, the 
field lines would be pulled by the accreting material into a split magnetic monopole, which is 
not present in any of the simulations discussed so far. Indeed, at small radii where the Hall 
spin- up is most efficient, the inward advection of poloidal magnetic field is almost exactly 
balanced by the outward magnetic diffusion enabled by the Hall effect. This balance is 
illustrated in the right panel of Fig. |31 which shows that, inside a radius ~ 10^^ cm, the 
radial component of the Hall-induced field line drift velocity on the equator has almost the 
same magnitude as the infall speed but with an opposite sign. The radial magnetic diffusion 
allows material to fall inward without dragging field lines along with it. The situation is 
similar in the azimuthal direction, where the Hall-induced field line drift speed and the fiuid 
rotation speed are nearly identical in magnitude but opposite in sign. The azimuthal Hall- 
induced drift prevents the field lines from being wound up continuously by rotation. The 
required toroidal current is provided by the bending of the field lines inside the disk (see left 
panel of Fig. 1; contrast with the unbent field lines in Fig. 7 of KLSIO). 

An interesting feature of the envelope collapse is that the infall speed of the equatorial 
material remains well below the free-fall speed. The equatorial collapse is slowed down three 
times (see right panel of Fig. [3]), for different reasons. The first slowdown near ~ 4 x 10^^ cm 
corresponds to the edge of the magnetic bubble infiated by field twisting (outside the region 
shown in the right panel of Fig. [3]), where a magnetic barrier forces the collapsing material 



over a large solid angle into a narrow equatorial channel (see Fig. 2 of iMellon fc Lil 12008 
and associated discussion). The second slowdown occurs near ~ 8 x 10^^ cm (or about 
500 AU), inside which the accreted magnetic fiux is left behind due to Hall diffusion. It 
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Fig. 3. — Left panel: map of the Hall-induced rotation speed in an initially non-rotating 
envelope at a representative time t = 2 x 10^^ s (Model NoROT; the unit for length is 
AU and for speed in the colorbar is cm/s). Also plotted are magnetic field lines (yellow), 
isodensity contours (black contours, in steps of 10), and velocity vectors (white). Right 
panel: infall and rotation speeds on the equator. Also plotted are the negative of the radial 
and toroidal component of the Hall-induced field line drift velocity, —i%r (dotted line) and 
~'f-%(f> (dashed). 
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is similar in origin to the ambipolar diffusion- induced accretion shoclc first discussed in 



Li fc McKed ( 1l996l : a similar behavior is also present in the Ohmic dissipation-only case, 
LKSll). The increase in the poloidal magnetic field strength enabled by Hall diffusion in 
the radial direction makes it easier to spin up the equatorial material through the Hall- 
induced magnetic torque. The third slowdown near ~ 10^^ cm is obviously due to centrifugal 
effect, since an RSD is formed interior to it. 

To check how robust the above results from Model NoROT are, we did a number of 
variants of the model. First, we reduced the initial field strength to 30% of its original value, 
from 35.4 /iG to 10.6 /iG (Model LoB). The result is shown in the left panel of Fig. H] (see 
dashed lines). An RSD is again formed at the representative time t = 2 x 10^^ s, despite the 
reduced initial field strength. It is somewhat smaller than that in the original case, which is 
to be expected since the Hall spin-up depends on current density, which increases with field 
strength. Another difference is that the infall speed outside the RSD is higher in the weaker 
field case, because the magnetic forces that oppose gravity are weaker. Second, we reduced 
the coefficient Q by a factor of 10, from 8.46 x 10^^ to 8.46 x 10^° in Gaussian cgs units 
(Model LoQ). An RSD is still formed at t = 2 x lO^^g ^j^i^ ^ radius of ~ 3 x 10^"^ cm or 
20 AU, it is much smaller than the RSD in the original Model NoROT. If we reduce the value 
of Q by another factor of 3 (to 2.82 x 10^° cm^ s"^ G"\ Model LoQl), an RSD still forms, 
but it extends barely outside the inner boundary of our computation domain, which has a 
radius of 10 AU. When the value of Q is reduced further to 8.46 x 10^^ cm^ s~^ G~^ (Model 
LoQ2), the RSD disappears. We therefore take 3 x 10^° cm^ s^^ G^^ as a rough estimate for 
the value of Q needed for the Hall-enabled formation of a relatively large RSD of tens of 
Alls or more in size. 

We have verified that the RSD formed in Model NoROT is little affected when the 



classical Ohmic resistivity used in KLSIO (based on iNakano et al.l |2002| ) is included. An 
enhanced resistivity can, however, reduce the size of the RSD significantly. Experimentation 
shows that it suppresses the Hall-enabled RSD formation in Model NoROT altogether when 
the resistivity is on the order of 10^°cm^s~^ or larger. The negative effect of the resistivity 
on the Hall spin-up is not surprising, since it reduces the electric current density in general, 
and the toroidal current density in particular; the latter is the driver of the spin-up through 
the Hall-induced magnetic torque. 

How does the critical value of Qc ~ 3 x 10^° cm^ s^^ G^^ estimated for RSD formation 
compare with the microscopic values of expected in realistic dense cores? In a separate 
study (LKSll), we have computed the values of Q as a function of hydrogen number density 
"^n = p/(2.33 X lO'^^g) , based on a simpl e prescription for the magnetic field strength 



as a function of density (INakano et al.ll2002l ) and a simplified chemical network for charge 
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Fig. 4. — Left: Comparison of infall and rotation speeds on the equator for the non- rotating 
Model NoROT (sohd hnes) and its two variants: Models LoB (dashed) and Model LoQ 
(dotted). Right: Comparison of the estimated critical value for the Hall coefficient Q^. for 
RSD formation (solid line) with the microscopic values of |(5| computed for the MRN grain 
size distribution (dotted) and 1 /xm sized grains (dashed), assuming a cosmic ray ionization 
rate of 10"^^ s'^ 
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densities ( iNishi et al.lll99l[ ). i ncluding two extrem e grain size distributions: (1) the standard 
MRN power-law distribution ( iMathis et al.lll977l : appropriate for diffuse interstellar clouds), 
and (2) grains of single, large size of 1 fim. These two cases should bracket the situation in 
dense molecular cloud cores, where some grain growth is expected. As can be seen from the 
right panel of Fig. HI although the microscopic values of Q can be larger than the critical 
value Qc at relatively low densities, they are smaller than Qc by about an order of magnitude 
at densities of order nn ~ 10^° cm~^ (or higher) that are crucial for RSD formation (see the 
left panel of Fig. [T]) . 



4. Conclusion and Discussion 



We have studied the collapse of a non- self-gravitating rotating, magnetized envelope 
onto a central stellar object of fixed mass to illustrate the influence of the Hall effect on 
disk formation. In this idealized problem, the formation of a rotationally supported disk 
(RSD) is completely suppressed by magnetic braking in the ideal MHD limit. Including 
a Hall term with a coefficient Q ^ 3 x lO^'' cm^ s~^ in the induction equation enables 
an RSD of radius > 10 AU to form. The RSD is formed not because the Hall effect has 
reduced the efficiency of the magnetic braking. Rather, it is produced because the Hall effect 
actively spins up the inner part of the equatorial material to Keplerian speed. The spin-up 
comes about because the collapsing envelope drags the magnetic field into a highly pinched 
configuration near the equator, producing a large toroidal electric current density, which 
forces the field lines to rotate differentially due to the Hall effect. The resulting twist of field 
lines yields the magnetic torque that spins up the equatorial material, even in the absence 
of any initial rotation. The spin-up is most effective in the inner part of the accretion flow 
where radial Hall diffusion enables the magnetic flux that would have been dragged into the 
central object by the accreted material in the ideal MHD limit to stay behind; the resulting 
pileup of magnetic flux at small radii is similar to the cases with either ambipolar diffusion 
or Ohmic dissipation. When the fleld direction is flipped, the Hall-induced magnetic torque 
changes direction as well. An implication is that the direction of the angular momentum of 
the material close to a protostar, including RSD, may depend more on the magnetic fleld 
orientation than on the initial rotation of the dense core that the star is formed out of, at 
least in principle. 

In practice, it is uncertain whether the Hall-induced magnetic torque can produce a 
sizable RSD or not. First, the value of the coefficient Q required to produce an RSD of 
tens of AUs in size in our idealized model appears to be larger than the microscopic values 
expected in dense cores, by roughly an order of magnitude. Whether the coefficient can be 
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enhanced somehow, perhaps through anomalous processes, is unclear. More importantly, the 
Hall effect is only one of the non-ideal MHD effects that are present in the lightly ionized, 
magnetized, dense cores. It can be dominated by ambipolar diffusion at low densities and 
by Ohmic dissipation at high densities. It will be interesting to explore the interplay of all 
these three non-ideal MHD effects and how they affect the collapse of dense cores and the 
formation of protostellar disks. Nevertheless, we have demonstrated that the Hall effect is 
unique among the non-ideal effects in its ability to actively spin up a magnetized collapsing 
flow, potentially to Keplerian speed, providing a new mechanism for disk formation. 

We acknowledge support by the Theoretical Institute for Advanced Research in Astro- 
physics (TIARA), by the National Science Council of Taiwan through grant NSC97-2112- 
M-001-018-MY3, and by NASA through NNG06GJ33G and NNX10AH30G. 
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